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ABSTRACT 

Linear theory provides a reasonable description of the velocity correlations of biased 
tracers both perpendicular and parallel to the line of separation, provided one accounts 
for the fact that the measurement is almost always made using pair-weighted statistics. 
This introduces an additional term which, for sufficiently biased tracers, may be large. 
Previous work suggesting that linear theory was grossly in error for the components 
parallel to the line of separation ignored this term. 
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- large scale structure of the universe 



> 

m 
o 

(N 

oo 
o 



1 INTRODUCTION 



Experiments will soon measure correlations in the peculiar 
velocity field traced by massive galaxy clusters. The linear 
theory description of these correlations as traced by dark 
matter particles whose motions are driven by gravitational 
instability from an initially Gaussian random field, has been 
available for about 20 years (Gorski 1988). The extension 
of this formalism to biased tracers, within the framework of 
perturbation theory, is in Fisher (1995); Sheth et al. (2001b) 
derive similar results within the context of the Halo model 
(see Cooray & Sheth 2002 for a review); and Regos & Szalay 
(1996) study this problem in the context of peaks (see Des- 
jacques 2008 for corrections to their expressions). All these 
analyses yield consistent results. 

However, measurements of the velocity correlations of 
massive haloes in simulations have given the impression 
that, although biased linear theory provides a reasonable 
description of the correlations between the velocity compo- 
nents that are perpendicular to the line of separation 
it is wildly discrepant for the components which are parallel 
to the line of separation (Croft & Efstathiou 1994; Peel 
2006), except on scales larger than about 100 Mpc. This 
discrepancy has been attributed to nonlinear effects, such 
as those described by Colberg et al. (2000) and Sheth & 
Diaferio (2001). 

The main purpose of this short note is to show that, in 
fact, linear theory does indeed provide a good description 
of provided one accounts for the fact that the measure- 
ment is actually pair weighted, and so the mean streaming 
motions may not be negligible. For dark matter, this mean 
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is sufficiently small that it can be ignored, but ignoring the 
mean streaming of massive halos towards one another is a 
very bad approximation. We show that keeping this term 
in the theoretical calculation provides substantially better 
agreement with the measurements. 

Section 2 provides the linear theory expressions for the 
velocity statistics of interest, Section 3 presents a compari- 
son of these expressions with simulations, and a final section 
summarizes our results. An Appendix provides some tech- 
nical details of the calculation. 



2 BIASED LINEAR THEORY 

In linear theory for dark matter, the mean approach velocity 
of particle pairs along the line of separation is given by 
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where H is the Hubble constant, /(fi) » Sl 5/9 for flat ACDM 
models, 
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and P(k) is the linear theory power spectrum of the density 
fluctuation field. This is often called the streaming motion; 
the term 1 + £ in the denominator reflects the fact that 
the average is over all pairs with separation r. For what is 
to follow, it will be convenient to write this as Wj™^) = 
(«j] m (r)), to emphasize the fact that the motion is parallel 
to the separation vector. In linear theory, the mean motion 
perpendicular to the line of separation is zero: (t>i m (>")) ~ 0- 
Dark matter halos are biased tracers of the underly- 
ing density field. If this bias is linear, but the velocities are 
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unbiased, then 



[l + b 2 ar)} <«n(r)> 



bvf?(r) [l + «r)], 



where & is the linear bias factor (see Appendix). 
The linear theory velocity dispersion is 
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where the smoothing window is related to halo mass by set- 
ting W{x) = (3/x)ji(x) with R = (3M/47rp) 1/3 , where p 
is the comoving number density. The factor of 3 makes this 
a one-dimensional velocity dispersion. The presence of W 
makes the velocity dispersion decrease slightly with increas- 
ing halo mass (Sheth & Diaferio 2001). Note that there is 
no additional dependence on halo bias factor. 
The two-point velocity correlation tensor is 



#ij(r) = (vi(x)vj(x + r)) 

= ^±(r)Sij -, + [*n(r) - *L(r)]fifj, 



(5) 



where r = |r|, r = r/r, and ^\\(r) and ^x(r) are the 
radial and transverse correlation functions (Monin & Ya- 
glom 1975). In linear theory, the velocity field is potential, 
so *||(r) = dr*x(r)/dr (e.g. Gorski 1988). For Gaussian 
initial conditions, the linear theory correlation between ve- 
locity components perpendicular to the line of separation 
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whereas the linear theory correlation between velocity com- 
ponents parallel to the line of separation is 
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A number of previous analyses have ignored the final term 
in this expression. The Appendix shows why, if the velocity 
correlations are estimated using pairs of tracer particles (as 
is commonly done), it should be included. 



3 COMPARISON WITH SIMULATIONS 

We compare the expressions above with measurements made 
in the Hubble Volume simulation (Evrard et al. 2002). The 
background cosmology for this simulation is flat ACDM with 
Qo = 0.3, a 8 = 0.9 and H = 100ft km s" 1 with ft = 0.7. We 
present results for halos within a very narrow mass bins 
centered on log 10 M jhT 1 Mq — 14, 14.5 and 15. There were 
28956, 21461, and 9352 haloes in each bin. 

The measured one-dimensional velocity dispersions of 
the halos are a v — 313, 303 and 283 km s" 1 , whereas the 
values predicted by equation (4) are slightly smaller: 289, 
276 and 258 km s" 1 . This may be an indication that non- 
linear effects of the sort discussed by Colberg et al. (2000) 
and Sheth & Diaferio (2001) have affected halo velocities. 
If we were to set W —* 1 in equation (4) as a crude way 
of accounting for nonlinear evolution, then this would make 
a v — 322 km s -1 . 

Figure 1 shows the mean streaming motions of massive 
halos along the line of separation. Notice that these velocities 
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Figure 1. Mean streaming motions of halos with masses 
log 10 (M//i- 1 M Q ) = 14,14.5 and 15 (bottom to top). Smooth 
curves show the linear bias prediction, bvi2(r), with b = 1.2,2 
and 3. 



can be large: the important point in the present context is 
that, on small scales, vi% S> (T v . 

The smooth curves show vf™(r) [1 + £(r)] of equa- 
tion (1), multiplied by bias factors of b = 1.2, 2 and 3. They 
provide a reasonably accurate description of the measure- 
ments. One might have thought that a more appropriate 
model would not multiply by 1 + £, and would include a fac- 
tor of W 2 (kR) in the integral of equation (1), for the same 
reason that it is included in the expressions for cr v , \&m and 
*x • We have found that including such smoothing produces 
too little streaming (compared to the simulations) on small 
scales. If the larger streaming on small scales is a conse- 
quence of nonlinear effects, then it may be that removing 
the smoothing term is a convenient way to approximately 
account for nonlinear evolution. (E.g., setting W — » 1 in 
equation 4 brings a v into slightly better agreement with the 
measurements.) Multiplying by 1 + £ serves a similar pur- 
pose. 

Figure 2 compares the predicted velocity correlations 
for these halos (curves) with the measurements (symbols). 
The solid curves show the predictions of equations (6) 
and (7). The dashed curves in the panels on the left show the 
result of ignoring the final term in equation (7); notice how 
the dashed curves grossly overestimate the measured signal 
(filled symbols) . This is the apparent failure of linear theory 
noted by Croft & Eftathiou (1994) and Peel (2006). How- 
ever, they are in good agreement with the measurements if, 
in the measurements, we add back the contribution from the 
mean streaming term (open symbols). This, and the reason- 
able agreement between the solid curves and solid symbols, 
indicates that linear theory is, in fact, not so badly flawed. 
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Figure 2. Velocity correlations of halos with mass log 10 (M/h~ 1 Mq) m 14, 14.5 and 15 (top to bottom) along (left) and perpendicular 
(right) to the line of separation. Filled circles show the measured correlations, open circles in the panels on the left show the result of 
adding the contribution of the mean streaming motions, {^i2(f*)) 2 /4, to the filled circles. Solid curves show the linear theory predictions; 
dashed curves in the panels on the left show the linear theory prediction when the mean streaming is assumed to be negligible. 



4 DISCUSSION 

The mean streaming motions of massive halos can be as 
large as 1000 km s _1 on scales of order 10/i -1 Mpc. In con- 
trast, the one dimensional velocity dispersions of these ha- 
los is of order 300 km s _1 . Nevertheless, these motions are 
rather well described by biased linear theory. The fact that 
the mean streaming motions of halos can be substantially 
larger than the halo velocity dispersions has an important 
consequence: they must not be ignored when making the 



linear theory estimate of velocity correlations ^ and &± 
for pair-weighted statistics. Ignoring this contribution leads 
one to predict that 'Pii is positive; i.e., objects stream along 
with one another, whatever their mass. The measurements 
show that while this is true for the lower mass halos (they 
are still quite massive!), is negative for the most mas- 
sive halos. This means that the most massive halos move 
towards, rather than along with, one another; this is consis- 
tent with linear theory provided one accounts for the fact 
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that the measurements are pair-weighted, so one should not 
ignore the contribution of the mean streaming motions. 

We note that, although they do not say so explicitly, 
halo model analyses of halo motions (Sheth et al. 2001a) 
have correctly included the effects of this mean streaming. 
It is the appearance of more recent studies which ignore 
this effect (e.g. we had communicated the importance of 
this term to Peel 2006, but his figures do not include it) 
which prompted our work. Because our analysis shows that 
linear theory can be used down to scales which are signif- 
icantly smaller than previously thought, we hope that our 
analysis will aid in the interpretation of data from the next 
generation of surveys. 

For deep surveys, the gain in accuracy is likely to be 
modest. This is because sky surveys typically estimate the 
correlation between radially projected velocities, separated 
by some angle 6 on the sky. This quantity is related to ours 
by 



#12 = #±cos# + (#n - #x)/(0,n,r 2 ), 



(8) 



where #| = *J- = a t V'-L, both are evaluated at scale 
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This can also be written as 
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#12 (r) = #x(r) Af- sin 9 + 

(10) 

showing that our extra term is most important for pairs 
which are close in angle but lie at different distances along 
the line of sight. For very deep surveys (ri and r 2 greater 
than 100 Mpc, say) the observable is dominated by #x for 
all scales where the addition of our extra term is important. 
Nevertheless, we emphasize that the term is present, and 
its inclusion greatly improves comparison with simulations, 
thus strengthening ones confidence in the validity of linear 
theory. 

Although accounting for mean streaming motions re- 
sults in substantially better agreement with the simula- 
tions, there remains room for improvement. Our treatment 
of how (f||(r)) should depend on halo mass on small scales 
is rather ad hoc. Methods motivated by perturbation the- 
ory are beginning to provide more detailed prescriptions for 
this term (e.g. Smith et al. 2008). Also, in perturbation the- 
ory, the power spectrum of peculiar velocities is expected to 
be slightly suppressed relative to the linear theory predic- 
tion, even on relatively large scales (e.g. Bernardeau et al. 
2002; Cooray & Sheth 2002; Pueblas & Scoccimarro 2008). 
Although including this effect is beyond the scope of this 
work, we note that accounting for it will lower the theory 
curves slightly, reducing the small discrepancy in the pan- 
els on the right of Figure 2. Incorporating these refinements 
into our analysis is the subject of work in progress. 
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APPENDIX A: LINEARLY BIASED TRACERS 

Dark matter halos are biased tracers of the underlying den- 
sity field. If this bias is linear, but the halo velocities are un- 
biased, then pair-weighted (as opposed to volume weighted) 
velocity statistics of halos may be biased relative to those of 
the mass, with the effect increasing with spatial bias. 

To see this, it is helpful write the expressions in the 
main text explicitly as averages over pairs. In linear theory, 
the mean streaming motions are 
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and it is understood that r = \ r i — r 2 \ . The factors of 1 + bi Si 
and l + b 2 5 2 represent the halo counts at positions n and r 2 , 
so their product is the pair-weight, and the averages are over 
all pairs separated by r. The algebra above follows because 



Velocity correlations of clusters 5 



(vi — v 2 ) = 0, and because, in linear theory, (SiVi) — 
and {6162V 1) — 0. On large scales the difference in the pair 
counts is small (£(r) <C 1), and V12 — ► fcfis^M where & = 
(61 + b 2 )/2. If the statistic is measured for a range of halo 
masses, then one simply sets both b\ and b 2 equal to the 
average bias factor b in the expressions above; this is the 
form we have used for the final expression. Note that V12 is 
biased relative to the dark matter, even though the velocities 
themselves are explicitly unbiased. 

Similarly, in linear theory, the mean correlation of halo 
line-of-separation velocities, averaged over halo pairs, is 

{(l + bi6i)(l + b26 2 )vi\\V2\\) 
<(1 + Mi)(l + M 2 )> 

{Vl\\ «2||) + (Ml b 2 S 2 fl|| V2\\) 
1 + &1&2 (6162) 

V2\\) + &iM<5i<52)(«i||«2||) + bib 2 {6iv 2 \\){5 2 v 1 n) 
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Note that, in the second and third lines, we have assumed 
Gaussian statistics to neglect three point averages, and to 
write the four point average as a product of two-point aver- 
ages. The final expression is what one obtains if a range of 
halo masses, with average bias factor 6, contributes to the 
statistic. On large scales, t>i™ and £(r) are both small, so 
"I'll w (fi|| v 2 \\). However, the second term in the expression 
may be non-negligible when b 3> 1; neglecting it will lead to 
unnecessary inaccuracies. A similar analysis of 9± does not 
yield an extra term because {6 2 vi±} = 0. 



